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Abstract. 

The possibility of unifying dark-matter and dark-energy has recently attracted 
considerable interest. In this so called quartessence scenario, a single component is 
responsible for both the clustering of matter and the accelerated expansion of the 
universe. A model archetype for such scenario is provided by the Chaplygin gas. 
Although this model is in agreement with the data on the expansion history, problems 
arise in the power spectrum of density fluctuations for adiabatic perturbations. In 
this contribution we consider other quartessence models and confirm that instabilities 
and oscillations in the matter power spectrum are a characteristic of more generic 
quartessence models, namely those with a convex equation of state. We show that, 
as in the Chaplygin case, this kind of problem can be solved by considering intrinsic 
non-adiabatic perturbations such that, as an initial condition, the perturbed fluid is 
gradient pressure free. We also discuss how the problems of adiabatic quartessence can 
be circumvented by other types of equations of state. 
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1. Introduction 

The twentieth century witnessed the establishment of a cosmological model consistent 
with a number of astronomical observations. In this standard model, the dynamics of 
the universe is governed by the local distribution of matter-energy, following Einstein's 
theory of gravitation. On large scales the universe is well described within a nearly 
flat Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmological scenario, and is 
presently undergoing an accelerated expansion. Despite the success of this model, several 
key features of the universe remain unclear. For example, what constitutes the dark- 
matter (DM) that would trigger gravitational clustering, and the dark-energy (DE) 
that is supposed to power the accelerated expansion? Are these two dark components 
separate entities or different manifestations of a single matter-content? 

Although there are some candidates for the DM from particle physics pQ, there is 
yet no evidence for these particles in laboratory experiments. For the DE the situation 
is even more puzzling since there is no natural candidate from fundamental physics, 
although the cosmological constant and a dynamical scalar field are the most popular 
[2j. Thus, the current information about DE and DM comes solely from astronomical 
observations. From a simplicity viewpoint it is interesting to investigate the possibility 
of describing the phenomenology associated with both DE and DM through a single 
matter-component, unifying dark matter (UDM) — or simply quartessence — [Sj, 
reducing to one the unknown components of the material substratum of the universe. 

A few candidates for such UDM appeared in the literature in recent years, the 
most popular being the generalized Chaplygin gas |11 , , IHl El E] , in what we call the 
Quartessence Chaplygin Model (QCM). For this specific model, several confrontations 
with observational data were performed and constraints were set on the model 
parameters. The generalized Chaplygin gas (as quartessence) appears to be compatible 
with all available data regarding the expansion history (see e.g. ref. [HI). However, 
for adiabatic perturbations, it was shown that the mass power spectrum presents 
strong oscillations and instabilities, allowing only models very close to the ACDM limit 
[TU] . This problem can be avoided if non-adiabatic perturbations with a specific initial 
condition {5p = 0) are considered jTTj. In [TU] it was argued that these oscillations and 
instabilities would be present in any quartessence model such that p = p{p). We argue 
that this is indeed the case for convex equations of state. This is illustrated through the 
investigation of two extreme models of convex quartessence, one with a very steep and 
another with a very gentle equation of state. We show that for these models, as in the 
QCM case, the above mentioned problems in the mass power spectrum can always be 
avoided with the same kind of entropy perturbation. On the other hand, we argue that 
the problems of adiabatic quartessence could be avoided by other types of equations of 
state. 

The paper is organized as follows. In section |21 we review the QCM and introduce 
the two mentioned extreme quartessence models. The treatment of linear perturbations 
in these models is discussed in section |3] The computation of the power spectrum 



Matter Power Spectrum for Convex Quartessence 



3 



in the adiabatic and nonadiabatic cases and comparison with observational results are 
developed in section |3] Finally, we sum up our results and present concluding remarks 
in sectional 



2. Phenomenological Models of Quartessence 

In this work we model the UDM component as a fluid with isotropic energy-momentum 
tensor, whose only independent components are the pressure p and the density p. We 
further assume that the space averaged dynamical variables are related by the equation 
of state (EOS) p = p{p), i.e. the background is a perfect fluid. For this EOS to be a 
quartessence candidate, it has to fulfill some conditions. For example, neglecting the 
contribution of baryons, in the FLRW model, for the universe to undergo acceleration, 
one must have at present (p + 3p) < 0. Thus we require p to be sufficiently negative 
when p is small. On the other hand, in the past, for the UDM component to behave 
like dark matter we must have p ^ p, such that the universe was decelerating and the 
density scaled as (where a is the scale factor). A remarkable property of quartessence 
is that the weak energy condition is never violated in such models. From the energy 
conservation equation 

p = -{p + p)3^ = -p{l + w (p)) 3^ . (1) 

it follows that the line p = —p cannot be crossed. Since at early times we should have 
w ~ 0, a straightforward prediction of these models is, therefore, that w > —1. If 
observations find that w; < — 1, this automatically rules out quartessence models in 
which the background is a prefect fluid (i.e. in which w is a function only of p) and 
where the dark sector is not coupled to other components. 

A simple example of an equation of state satisfying the above conditions is given 
by an inverse power law jlj , O , jHl 13 E] , 

7\^4{a+l) 

p = — — (Chaplygin Quartessence), (2) 

where M has dimension of mass and a is a dimensionless parameter. The energy 
conservation equation (0) has a simple analytic solution for this EOS 

3(a+i) n 



PCh — PChO 



(3) 



Here Oq is the present value of the scale factor and A = (M^/pcho)^"^^''. As expected, 
when a/ao ^ 1, we have pch oc a"^ and the fluid behaves as CDM. For late times, 
a/flo ^ 1) and we get pch = —Pch = —M^ = const, as in the cosmological constant 
case. 

The QCM has been extensively discussed in the literature recently, both setting 
observational limits to the model and studying its properties and physical motivation. 
Nevertheless, since the QCM is a prototypical and particularly simple example of 
quartessence we shall briefly review a few results concerning the parameters of the 
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model. This will help to fix the notation and gain some intuition, as several aspects are 
similar in generic models of quartessence. 

Since in QCM the equation of state parameter w is a monotonic decreasing function, 
it is easy to see from equation (jT)) that the equation of state goes asymptotically to a 
cosmological constant form p (p) = —p. The minimum value of the density, obtained 
from the solution of this equality, will be denoted by pmin and is simply given by in 
the QCM. Naturally, the current density pcho cannot be lower than pmin (otherwise, as 
the p = — p line cannot be crossed, QCM would always have had w < —1, never acting 
as dark-matter). Thus pmin < Pcho < oo, which restricts the parameter A to the range 
< a < 1. Clearly for A = the fiuid behaves as dust, while for A = 1 it acts as a 
cosmological constant. 

Now, let us discuss the range of allowed values for a. In most of the analyses in 
the literature it is assumed that < a < 1. These limits are imposed because adiabatic 
perturbations of the QCM are implicitly assumed. Notice that, for these models, there 
is an absolute maximum value for the adiabatic sound speed = [dp/ dp), which occurs 
for p = Pmin and is given by c^^ax — ^- Therefore, to avoid superluminal propagation 
of signals the upper limit (a < 1) is imposed in the adiabatic case. Also, demanding 
stability of adiabatic perturbations implies {dp/ dp) > 0, which requires a > 0. However, 
as shown by Sandvik et al. fU]; adiabatic QCM is ruled out by current data unless 
a is very close to zero (corresponding to the ACDM limit). Nevertheless, if entropy 
perturbations such that the effective sound speed vanishes are allowed, the model is 
again consistent with observational data In this nonadiabatic case, the above 

limits on a are unnecessary. For instance, if a > 1 there is no problem of superluminal 
propagation for these specific entropy perturbations and even models with a < can, 
in principle, be in accordance with observations. Since for quartessence we should have 
u; ^ as p ^ oo, values of a < — 1 are not allowed. Therefore, here we only require 
that a > —1. 

The QCM power spectrum exhibits strong oscillations or instabilities, when 
adiabatic perturbations are considered and a is not very close to zero. As we shall see in 
the next section, this is due to the development in the fiuid of a non null adiabatic sound 
speed at recent times. For QCM the adiabatic sound speed is related to the equation 
of state parameter by = a{M'^ / p)^"~^^^ = —aw. To avoid such oscillations in linear 
scales, the parameter a of the QCM has to be restricted to |a| < 10^^, very close to the 
ACDM limit POj- One may wonder if a steeper variation of the equation of state would 
reduce the oscillations, allowing a wider range of parameters, distinct from the ACDM 
model. Let us consider for example an exponential equation of state 



In this case, the adiabatic sound speed has an additional exponential suppression, 
cj. = a exp {—ap/M^). Thus, one could expect that high ratios of p/M^ would attenuate 
the oscillations in the power spectrum, without imposing strong constraints on the 
parameter a. However, the equation of state parameter w also has the same additional 




(Exponential Quartessence) . 



(4) 
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exponential factor and, if the ratio p/M^ is always high, the Universe does not enter in 
an accelerated phase. Therefore, as it will be shown explicitly in the next section, this 
modification in the equation of state does not help in solving the problem and indicates 
that it is a general property of the model. 

Before introducing our third type of quartessence, let us discuss some noteworthy 
values of the exponential model parameters. First, for Eq. to be a quartessence 
candidate (i/; ~ as p ^ oo) the parameter a must be non negative. The condition 
a > also guarantees stability for adiabatic fluctuations. The minimum energy density 
of this model is given by pmm = M'^W {a) /a, where W {x) is the Lambert function 
defined by W (x) exp {W {x)) = x ^2]- Therefore, the maximum adiabatic sound speed 
is cLax = W{a). 

Another extreme case, such that the model has a ACDM limit, is given by: 

p = — -p — (^-/M^)]"^ (Logarithmic Quartessence) . (5) 

In this case, for w < to be an increasing function of p, either a > or a < — e. 
The minimum density is given by pmm = exp (aVT and cj.^^^ = 1/W {a~^). 

As in the exponential and QCM cases, for a = we recover the ACDM model for the 
background. 

In the adiabatic case, where the condition a > is fulfilled such that > 0, the 
three models discussed above have a common property, namely that the equation of 
state is convex, i.e. dPp/dp^ = dcj./dp < 0. Therefore, the maximum value of the sound 
speed is reached at the minimum value of the density. We shall refer to models that 
share this property as convex quartessence. The above new Anzdtze (Eqs. |3] and El) 
cover two extreme cases of convex quartessence equations of state, one with a very steep 
and other with a very gentle variation of the pressure with the density. We believe that 
by analyzing these two examples we may have indications of generic behaviors for any 
(monotonic, well behaved, smooth, etc.) convex quartessence EOS. In the next sections 
we shall focus on the matter power spectrum of these models, show that they have the 
same oscillations as the Chaplygin case and check that the introduction of nonadiabatic 
fluctuations solves this problem. 



3. Evolution of Linear Perturbations 



The relativistic equations that govern the linear evolution of scalar perturbations in a 
multi-component fluid, in the synchronous gauge, are [ini Ej 



a' . /, h'r\ „ a' 



5[ + 3(4 - Wi)-5, = -(1 + w,)[kvi + ^\- ?,Wi-Tu (6) 



a' cl- . „ Wi 



„: + (l-34)-.. =^H. + ^ff., (7) 
h"r + -h'r = - V(l + ScDSnGpia^Si - 2A7rGa^y" PiT,, (8) 
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where 5^, Vi, and Fj are, respectively, the density contrast, the velocity perturbation, 
and the entropy perturbation of component i, is the trace of the metric perturbation, 
and k is the commoving wave number. For the sake of simplicity, we assumed that both 
the spatial curvature and the anisotropic pressure perturbation vanish and that the 
energy-momentum tensor of each component is separately conserved. In the equations 
above the prime denotes derivative with respect to conformal time {dt = a drj) and, as 
usual, c^j = p'i/p'i and Wi = pi/ pi are, respectively, the adiabatic sound speed and the 
equation of state parameter of component i. 

We will consider only two components in the universe: baryons {pi, = 0) and 
a quartessence component. To glimpse general features of quartessence, we consider 
the two extreme models introduced above: exponential (Eq. and the logarithmic 
quartessence (Eq. E}. As remarked, these models have the ACDM model as a limiting 
case for the background when a = 0. At first and higher orders in perturbation theory, 
however, these quartessence models with a = 0, have distinct behavior as compared to 
that of dark-matter in ACDM fT^ITT)] . 

The independent conservation of energy-momentum for baryons and quartessence 
leads to 

pb = pboa~^- (9) 

For the exponential and logarithmic cases, the energy conservation equation does not 

have an analytical solution. We rewrite Eq. as 

du u , . , , 

— + 3-(l + w)=0, (10) 
da a 

were u = p/ po, w = w (^u, a, , and A = M^/pQ. For the exponential case one has 
A f au\ 

w = exp — — (11) 



u V ^ . 
whereas the logarithmic case leads to 



«; = _^An^) . (12) 



u \ A, 

Therefore, given a and A we numerically solve (fTUI) with the initial condition m(1) = 1 
to get u{a) and then p = pqu. Notice that, for each a, A has to be chosen such that the 
condition < Pmin/Po < 1 is fulfilled. 

Let us write Eqs. (0111) for the baryon-quartessence system. For baryons, as c^^ = 
and Wb = 0, Eq. ((Zj) implies Vb oc a~^. As a matter of simplification, we take Vb = 0. 
Now Eq ^ implies h'j^/2 = —5^. Using this result in Eq. (jHJ and changing the "time" 
variable, we obtain, for the flat case {Qbo + ^qo = 1); 

(M) § - # ^ - [(^ ^ + ^'"'rj } .(13) 

For the quartessence component, Eq. ^ can be rewritten as 

— ^ + 3 (c -Wg) — = -(1 + Wq) { ^ - — ] - 3— Fg, (14) 
da ^ ^ a ^ \aa da J a ^ 
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whereas Eq. ((Tj) gives 



The derivatives of the scale factor with respect to t are given by 



a = i/oy^ + (l-^feo) (16) 

a = + (1 - fi.o) u{l + 3w,) . (17) 

In the next section we apply these results to compute the matter power spectrum for 
quartessence models. 



4. The Power Spectrum 

In order to obtain power spectra for the considered models we evolved Eqs. p3H15p 
from z = 500 to 2; = 0, using the solutions of Eq. (fTn|) . At z = 500 the quartessence 
fluid behaves like CDM and, to take this into account, proper initial conditions were 
established: we assumed Vq = 0, a scale invariant primordial spectrum, and used the 
BBKS transfer function ^Tj, with the following effective shape parameter [TH} I19j. 



eff = i^bo + ^U^)h exp -Qbo - —jj 



where the term VL^J^^ = (1 — fi^o) lima_>o(Ma^) is the effective matter density. 

In Fig. ^ we present the baryon and the total mass power spectra, for 
the exponential quartessence model, in the adiabatic case {Tq = 0), for a = 
0, 10-^10-^10-3 (from top to bottom), h = 0.7, ilbo = 0.04, and we choose M^/po 
such that all models have Fe// = 0.18. The total power spectrum can be written as 
P{k) oc IQbSb+^qSql"^ ■ The normalization of the mass power spectra is arbitrarily fixed at 
k = 0.01 h Mpc~^. The (blue) squares are data points from the Sloan Digital Sky Survey 
(SDSS) as compiled in [201, whereas the (red) circles correspond to the Two Degree 
Field Galaxy Redshift Survey (2dFGRS) from The presence of oscillations in the 
total power spectra, due to a non-vanishing adiabatic sound speed of the quartessence 
component, is apparent in the figure. Further, even for the baryonic component, the 
spectrum is compatible with the data only for a <^ 10~^. In Fig. |21 we display 
the baryon and total mass power spectra for logarithmic quartessence, for the same 
parameter values as in Fig. ^ Again, oscillations and power suppression are present in 
the total power spectrum. However, for a < 10~^, the baryonic power spectrum is still 
compatible with the data. 

Notice that the amplitude of the initial power spectrum could be reescaled such 
that the baryon spectrum provides a better fit to the galaxy data. However, the 
normalization of the total power spectrum is constrained by gravitational lensing 
observations. Therefore, as remarked in JO], it is not possible to reconcile the constraints 
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Figure 1. Baryon (left panel) and total (right panel) mass power spectra for 
exponential quartessence, in the adiabatic case. The curves from top to bottom 
correspond to a = 0, 10^^, 10^"*, 10~^ and Af'/po is chosen such that all models have 
T^ff = 0.18. The (blue) squares are the data points from the Sloan Digital Sky Survey 
as compiled in [50] and the (red) circles are from the 2dF galaxy redshift survey as 
compiled in PT] , 
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Figure 2. Baryon (left panel) and total (right panel) mass power spectra for 
logarithmic quartessence, in the adiabatic case. The curves from top to bottom 
correspond to a = 0, 10^^, 10^"*, 10^"^ and jp^ is chosen such that all models have 
Fe// = 0.18. The (blue) squares are the data points from the Sloan Digital Sky Survey 
as compiled in [201 and the (red) circles are from the 2dF galaxy redshift survey as 
compiled in pT] , 
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Figure 3. Baryon mass power spectra for the model with exponential quartessence 
(left panel) and logarithmic quartessence (right panel), in the non-adiabatic case, 
a = 0,0.1,0.2 and M'^/po chosen such that all models have Fe// = 0.18. The (blue) 
squares are the data points from the Sloan Digital Sky Survey as compiled in |20| and 
the (red) circles are from the 2dF galaxy redshift survey as compiled in |21 |. 



from both the galaxy power spectrum and gravitational lensing observations for non 
negligible values of a and adiabatic perturbations. 

As discussed in jTT] the oscillations (for > 0) and instabilities (for < 0) in 
the quartessence mass power spectrum that are present in the adiabatic case, have their 
origin in a non- vanishing value of the right-hand side in Eq. (jZj). To avoid these effects 
we choose, for the quartessence component 



Thus the condition (jl9|) is equivalent to set 6pg = 0. Moreover, it is sufficient that 
relation (fT^ be satisfied at some initial time t, since this condition is preserved along 
the linear evolution. This can be easily seen by differentiating ()20j) with respect to time 
and substituting (fTT^ . which implies that d6pq/dt = 0, with no consideration about the 
exact form of the fiuid equation of state. 

In Fig. ini we show the baryon spectrum for exponential and logarithmic 
quartessence, in the case = ~c^^6q/wg, for the same parameter values of Figs. [T] 
and 121 The results are visually indistinguishable from the predictions of the ACDM 
model. The same holds for the quartessence and, therefore, the total power spectrum. 

It is important to remark that the above considered entropy perturbations are not 
isocurvature ones, since there are initial density perturbations, although no pressure 
fiuctuations. Moreover, if one were to consider the evolution of the coupled baryon- 
photon fiuid, we would set Spb/pb = Spq/ Pq = {'i/^)5pr/ Pr, which is sometimes referred 




(19) 



^From the definition of entropy perturbation ^^-i we have 
5p = pT + 0^5 p. 



(20) 
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in the literature as "adiabatic perturbations". Therefore, the entropy perturbations 
needed to ehminate the instabihties and oscillations in the power spectrum are 
not straightforwardly ruled out by cosmic microwave background radiation (CMB) 
observations, which have discarded purely isocurvature modes. A complete study of 
the CMB spectrum in quartessence models with entropy perturbations such as given by 
Eq. (fT^ remains to be done and would be an important test for the model. 

5. Discussion 

In what is currently considered the standard cosmological model, two exotic components, 
DM and DE, are invoked to explain two different phenomena: clustering of matter and 
cosmic acceleration. In spite of their success in explaining most of the cosmological 
observations, the exact nature of these two components remains a mystery. In fact, so 
far, it has not been proven that dark matter and dark energy are two distinct substances. 
This unsatisfactory situation stimulated the search for viable quartessence (unified-dark- 
matter) models [^1^. 

In this work we argue that adiabatic quartessence models with a convex equation of 
state always suffer from oscillations and suppression in the power spectrum, rendering 
those models inconsistent with the data. This can be understood in very simple terms, 
since, for that class of model, the maximum sound speed always occurs at the minimum 
density (at which p = —p). Therefore, the epoch of accelerated expansion is also a 
period of high adiabatic sound speed. This situation is different for equations of state 
that change concavity. In this case, the moment of maximum adiabatic sound speed can 
be disconnected from the start up of accelerated expansion. Such models are currently 
under investigation and may provide a power spectrum in agreement with the data. 

Another alternative is to consider intrinsic entropy perturbations. In Ref. ^T] it 
was shown that, for QCM, if Tg = —c1g6g/wq, or equivalently, that initial conditions are 
chosen such that 6pq = 0, the oscillations and instabilities present in the adiabatic case, 
disappear. In this work, we have shown that the same type of entropy perturbations 
also solves the problem of the power spectrum for both the exponential and logarithmic 
quartessence models. This occurs because, in this case the effective sound speed vanishes 
(^seff — ^p/^P = 0), such that there is no gradient term in the "Euler equation" ((7j), and 
thus no oscillations and divergences appear in the power spectrum. We therefore argue 
that such entropy perturbations render generic quartessence models in agreement with 
power spectrum observations. However, although the nonadiabatic quartessence power 
spectrum is indistinguishable from that of "concordance" models like ACDM, their 
behavior beyond the linear regime is quite different and, in principle, can be clearly 
distinguished with weak lensing observations [TK] . 

To conclude, we remark that the quartessence scenario offers an alternative to 
concordance models, providing distinct predictions that can be tested with current data. 
Several dark energy models, being as "exotic" as quartessence, are nevertheless much 
more degenerate with the ACDM model. It is therefore desirable to put some effort into 
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deriving observational consequences of quartessence models. 

Adiabatic convex quartessence now seem to be ruled out by the data. Nonadiabatic 
convex models might be discarded by weak lensing. However, models in which the 
equation of state changes its concavity are a promising alternative, even for adiabatic 
perturbations [SB]- Such models deserve to be further studied, offering potential 
alternatives for dark-matter and dark-energy unification. The search for quartessence 
models provides also a test for the robustness of the dark-matter plus dark-energy 
paradigm. If generic unifying-dark- matter models can be ruled out, this would give 
a strong support for the concordance cosmological models. 
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